Mon. Not. R. Astron. Soc. 000,[T]{8](2002) Printed 2 March 2013 (MN KTeX style file v2.2) 



On how leakage can affect the Star Formation Rate estimation using 
Ka luminosity 

M. Relano^*, R. C. Jr. Kennicutt^, J. J. Eldridge^ ^ J. C. Lee^ S. Verley^ 

^ Departamento de Fisica Tedrica y del Cosmos, Universidad de Granada, Campus Fuentenueva, Granada, Spain. 
^ Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CBS OHA. 
^ Department of Physics, University of Auckland, Private Bag 92019, Auckland, New Zealand. 
^ Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA. 



ABSTRACT 

We present observational evidence that leakage of ionising photons from star-forming regions 
can affect the quantification of the star formation rate (SFR) in galaxies. This effect could par- 
tially explain the differences between the SFR estimates using the far ultraviolet (FUV) and 
the Ha emission. We find that leakage could decrease the SFR(Ha)/SFR(FUV) ratio by up 
to a 25 per cent. The evidence is based on the observation that the SFR(Ha)/SFR(FUV) ratio 
is lower for objects showing a shell Ha structure than for regions exhibiting a much more 
compact morphology. The study has been performed on three object samples: low luminosity 
dwarf galaxies from the Local Volume Legacy survey and star-forming regions in the Large 
Magellanic Cloud and the nearby Local Group galaxy M33. For the three samples we find 
differences (1.1 — lAcj) between the SFR(Ha)/SFR(FUV) for compact and shell objects. Al- 
though leakage cannot entirely explain the observed trend of SFR(Ha)/SFR(FUV) ratios for 
systems with low SFR, we show the mechanism can lead to different SFR estimates when us- 
ing Ha and FUV luminosities. Therefore, further study is needed to constrain the contribution 
of leakage to the low SFR(Ha)/SFR(FUV) ratios observed in dwarf galaxies and its impact 
on the Ha flux as a SFR indicator in such objects. 
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1 INTRODUCTION 

Measuring the star formation rate (SFR) in galaxies has long been 
a major field of study in astronomy. Two of the most widely indica- 
tors of the star formation are the Ha recombination emission line 
and the UV non-ionising continuum. The Ha emission is produced 
by the recombination of the hydrogen ionised by the most mas- 
sive (M>20 Mo) and short-lived (t< 1-lOMyr) stars. The UV 
emission primarily originates in the photospheres of lower mass 
(M > 3 Mo) and long-lived stars (t < 200 Myr). Despite the differ- 
ent average timescales in the stellar lifetimes both indicators have 
been probed to trace well the SFR in galaxies over a range of lumi- 
nosities and environments ( Kennicutt|1998| ). 

The robustness of both indicators has been questioned for ob- 
jects whose star formation history consists of a set of starburst 
episodes rather than a constant SF R (^Sullivan et al.|2004 Iglesias-| 



paramo et al.||2004[ [BoselH et al.||2009[ |Bell & Kennicutt||2001| ). 
Recently, |Lee et al!] p009^ showed that, in a complete sample of 
^300 star-forming galaxies within 1 1 Mpc volume, extinction cor- 
rected Ha and far ultraviolet (FUV) emission give consistent SFRs 
for normal spiral galaxies (SFR ^ 1 Mo yr ~^). However, for low 
luminosity dwarf galaxies (SFR ^ 0.1 Mo yr~^) Ha emission 
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tends to under predict the SFR relative to the prediction of FUV 
emission. For galaxies with SFR < 0.003 Mo yr~^ the ratio of 
SFR(Ha)/SFR(FUV) is on average lower than expected by a factor 
of two. 

Several explanations have been suggested for explaining this 
discrepancy (see e.g. Meurer et al.|2009[|Lee et al.|2009| for a re- 
view): the amount of dust within star-forming regions, star forma- 
tion history (SFH), porosity of the interstellar medium, stochastic- 
ity of the Initial Mass Function (IMF) and IMF variations. 

The differences in the SFR derived from Ha and FUV for low 
star-forming systems are present even when Galactic foreground 
extinction corrections have been applied to the observed luminosi- 
ties. For normal spiral galaxies (see also |Botticella et al.|2011| ) the 
SFR derived from FUV tend to be lower than the one derived using 
Ha emission, but after accounting for internal attenuation agree- 
ment is found using both tracers. However, for low star-forming 
systems internal attenuation increases the discrepancies between 
the two SFR estimates ( [Meurer et al.|2009[|Lee'et al.|2009| ). 

A departure in the assumption of constant star formation im- 
plied in the calibration relations can affect the value of SFR de- 
rived from FUV or Ha emission: for a galaxy with systematic star 
formation bursts the deficiency of ionising short-lived stars rela- 
tive to the longer-lived lower mass stars emitting in the FUV will 
produce an Ha-to-FUV ratio lower than expected from a constant 
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SFH ( [Sullivan et aL|2004[[lgIesias-Paramo et al.|2004| ). However, 
to explain the observed variations in the SFR ratio the low star- 
forming systems must have very intense (a factor of 100) and very 
long (^lOOMyr) bursts of star formation (Iglesi as-Paramo et al. 
|2004| ). Such bursts of star formation do not seem to occur in these 
systems based on SFHs reconstructed from resolved stellar popula- 
tions ( |Weisz et al.|2008l|Mcquinn et al.|2009] ). 

Stochastic effects of an invariant IMF (e.g Cervino & Lurid^ 
|iana|2004l [Corbelli et al.|2009j or variation of the maximum stel- 
lar mass that can be formed in a stellar cluster leading to an inte- 
grated galactic initial mass function (IGIMF) (Weidner & Kroupa] 
[2005 , Pflamm-Altenburg et al. 2009 ), have been proposed to ex- 
plain the differences in the observed SFR(Ha)/SFR(FUV) ratios. 
The assumption of an IGIMF describing the stellar population of 
the galaxy seems to produce results that explain the observations 
reasonably well, but it underestimates the Ha luminosity in sys- 
tems with low SFR activity (Fumagalli et aL[[20lT| [Weisz et"ar] 
[2012| ). Besides, the variation of an upper mass limit of the IMF is 
still under discussion, as there is some evidence showing that this 
variation does not exist ( [Calzetti et al.|2010| ). [Eldridge] ( [2011[ ) and 
([Fumagalli et al.|20Tl] ) showed that stochastic IMF sampling com- 
bined with a cluster mass function and a cluster age distribution 
can well explain the differences in the SFR(Ha)/SFR(FUV) ratios. 
Using models of different SFHs and a fully populated IMF, [Weisz| 
|etari ( [20T2l ) compared the predictions of SFR(Ha) and SFR(FUV) 
with the observations and were also able to reproduce the decline 
of the observed SFR(Ha)/SFR(FUV) for lower luminosity systems. 
However, using the same models, [Weisz et aL[ ( [2012[ ) were not able 
to find a strong agreement between the Ha equivalent width or the 
UV luminosity modelled distributions and the observed ones. 

One of the other plausible explanations for the observed SFR 
ratio is the leakage of Lyman continuum photons from star-forming 
regions. In the conversion of Ha luminosity to star formation, it is 
assumed that all ionising radiation is absorbed by the neutral hy- 
drogen in the star-forming regions and therefore every Lyman con- 
tinuum photon will produce an Ha photon. However, it is known 
that some of the ionising photons can escape the regions of star for- 
mation (e.g. | Oey & Kennicutt| 1997 [[Relano et al. [20021 [Eldridge&l 
[Relano |201l| ) and ionise the diffuse interstellar gas (DIG) in nor- 
mal galaxies ( [Zurita et al. ,2000, |2002t , Wood et al. 2010) . [Oey eTaT] 
( [2007[ ) have shown that the fraction of diffuse ionised gas seems 
to be higher in dwarf galaxies than that observed in normal galax- 
ies. However, whether this translates into a difference in the escape 
fraction of ionising photons for each type of galaxy remains uncer- 
tain ( [Hanish et al.|2 010 ). If all the ionising photons escaping from 
the star-forming regions are able to ionise the DIG they should be 
included when integrating the total Ha emission from the whole 
galaxy. Therefore, for a constant SFH, we should not expect dis- 
crepancies between the SFR estimates from the total Ha or FUV 
luminosity of galaxies. However, if there are ionising photons es- 
caping the galaxy disk or those ionising the DIG are not detected 
( [Melena et ar][2009[ [Hunter et aT][2010[ l, then we expect to find 
discrepancies between the SFRs obtained from Ha and FUV lumi- 
nosities. Recently, Eldridge ('2011) has shown that leakage of ion- 
ising photons in combination with pure stochastic IMF sampling 
better reproduce the observed SFR(Ha)/SFR(FUV) ratio at very 
low luminosities. 

There is some evidence that the Hll regions with Ha shell 
morphology (normally called bubbles and superbubbles) have in 
general a higher escape fraction of ionising photons ( |Oey & Ken-[ 
[nicutt|1997] ), based on the idea that once the shells are formed the 
ISM can fragment and produces voids of gas where the ionising ra- 



diation can escape ( [Dove et al.|2000[ ). [Hoopes & Walterbos[ ( [2000[ ) 
tried to apply this idea to a sample of H II regions in M 33 and 
found no correlation of compact or diffuse star-forming regions 
with the expected leakage fraction of ionising photons. Recently, 
[Grossi etaTjpOlOl ) found that the star-forming regions in M 33 as- 
sociated with younger and lower-mass stellar clusters seem to have 
higher escape fractions of ionising photons. 

In view of the renewed interest in the consistency of Ha and 
UV-based SFRs we have undertaken a more comprehensive exam- 
ination of the possible role of ionising photon leakage from H II re- 
gions. We present here observational evidence that even if the frac- 
tion of ionising photons escaping from galaxies and star-forming 
regions is not able to completely explain the observational trends 
by itself, it should be taken into account when trying to explain 
the observed phenomenology in combination with other mecha- 
nisms. We estimate here that leakage could decrease significantly 
the SFR(Ha)/SFR(FUV) ratio. Therefore, in order to explain the 
differences in the quantification of the SFR based on FUV and Ha 
luminosities, one should constrain the magnitude of leakage from 
star-forming regions. 

Our study is based on a morphological classification of the 
Hll regions in low surface brightness dwarf galaxies in the Lo- 
cal Volume Legacy (LVL) survey sample ( [Dale et aL[[2009[ [Lee[ 
[et al.|2011[ ). Using the Ha images from the LVL archive, we per- 
form the morphological classification of the galaxies based on the 
morphology of the most luminous H II regions in each galaxy. The 
same morphological discrimination is applied on Ha images of in- 
dividual regions in the Large Magellanic Cloud (LMC) and the 
nearby galaxy M33. The most luminous Hll regions show Ha 
luminosities and SFRs comparable to the low surface brightness 
dwarf galaxies. Besides, the star formation for these galaxies occur 
in the form of short-lived burst episodes similar to those in the H II 
regions (see e.g. [Gerola et al.|1980| [Weisz et al.|2012[ ), therefore 
it is reasonable to study and compare both object samples. In Sec- 
tion 2 we explain the data for each sample and the morphological 
classification we have applied to the objects. In Section 3 we esti- 
mate the SFR(Ha)/SFR(FUV) for each object in the samples and 
show the trends with the Ha morphology. We analyse the results 
and show the importance of the leakage effect in the observational 
trends in Section 4 and the conclusions are presented in Section 5. 



2 METHODOLOGY 

2.1 Local Volume Legacy galaxy sample 

We classified the LVL galaxies from the Data Release f]( |Lee et"aL 
[2009[ ) in terms of morphology using the Ha images available in the 
archive. We chose only galaxies with RC3 type T 10 and 1 1 to avoid 
grand design spirals. The classification was done visually by one of 
the authors (MR) in terms of morphology of the H II regions: com- 
pact are galaxies presenting compact knots (one or several), mixed 
are galaxies presenting compact knots and filamentary structures, 
and shells are galaxies with clear shells (one or several). Ha im- 
ages of some examples of each classification are shown in the top 
row of Fig.[T] while the FUV images of the same objects are shown 
in the top row of Fig.|2] The integrated Ha and FUV fluxes for the 
LVL galaxy sample were presented in [Kennicutt et al.[p008j and 
in [Lee et ar[ ( [2010[ ), respectively. 
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The comparison of the Ha and FUV SFRs for the galaxy sam- 
ple was done in |Lee et aL] ( |2009| ) showing that for low luminosity 
dwarf galaxies Ha luminosity tends to underpredict the SFR rela- 
tive to the FUV luminosity. The trend is seen not only for Galactic 
extinction corrected SFRs but also when internal extinctions were 
applied. The Ha extinction was derived using the B aimer decre- 
ment while the extinction in the FUV was obtained using the total 
infrared (TIR) to FUV flux ratio and applying models from Buat 
|et al.| ([2005 ). For those objects with no available TIR emission es- 
timates scaling relations were used to derive the FUV extinction. 
Besides, scaling relations were used when Balmer decrements or 
TIR/FUV ratios were not available (see |Lee et al.|2009] ). 

2.2 Sample of star-forming regions in M 33 

We selected a sample of star-forming regions in the nearby galaxy 
M33. This galaxy is a particularly suitable object for this study 
because its distance (840 kpc; |Freedman et al.|1991| ) allowed us to 
observe the H II regions as individual objects and to perform a mor- 
phological classification of the sample. Besides, the large number 
of star-forming regions in the disk of M33 allowed us to create a 
significant statistical sample of objects. Some clear shells have been 
observed in several wavelength ranges from the optical to the Far 
Infrared (FIR) with the recent observations from Herschel ( |Kramer| 
[et al. 2010, Verley et al. 2010 ). An example is shown in the lower- 
right corner of Fig. [T] this shell is also clearly revealed in the FIR 
bands from Herschel observatio ns ([Verley et al.|2()T0| and will be 
studied in detail in |Relano et al.| ( [2012| in prep.). 

We used the continuum- subtracted Ha image of M 33 from 
[Hoopes & Walterbosl { [2000| ) to select our sample. The angular res- 
olution of this image, 5 corresponding to a linear scale of ^20 pc, 
is good enough to identify the H II region sample. The H II region 
selection was done visually by one of the authors (MR), as the aim 
was to find a set of relatively isolated objects with a clear mor- 
phology which facilitates the classification (see also Relano et al.| 
|2012| in prep.). The morphological classification was done using 
the wide-field KPNO mosaic Ha images from the Local Group 
Galaxy Survey (LGGS) collaboration (Mass ey et al.|2006[ ). These 
Ha images provide a much better spatial resolution (^0.8 arcsec, 
corresponding to ^4 pc) allowing a more accurate morphological 
classification. Our final sample consists of 117 objects. We applied 
the same criteria for the morphological classification as described 
in section [2?7] and found 10 compact regions, 46 mixed regions and 
61 regions with shell morphology. 

We performed photometry in the FUV, Ha, 24 /xm and 8 /xm 
bands using individual apertures for each H II region. The FUV im- 
age w as taken from the GALEX fMartin et al.|2005[ |de Paz et al.| 
[2007| data archive. The angular resolution for this image was ~4'^ 
The 8 /im and 24 /xm emissions were obtained from the IRAC and 
MIPS (Spitzer) images with ^^3^^ and 6" resolution, respectively. 
The data reduction of the Spitzer images is explained in| | Verley et al.| 
([2007 ). As the images were used to obtain the TIR luminosity of 
the objects, we needed to eliminate the stellar contribution in these 
bands. We used the IRAC-3.6 /xm image and followed the relation 
proposed by Helou et al. (2004) for this purpose. Finally, all the im- 
ages were smoothed to a common 6" resolution (the 24 /xm band 
resolution, which is the lowest one for the images of M 33 used in 
this study) and registered to the same pixel size as the Ha image 
from |Hoopes & Walterbosl ( [2000| . 

The Ha photometry was performed on the Hoopes & Walter- 
bos image, as the Ha image of the LGGS collaboration presented 
some saturated zones in the brightest part of the most luminous 



Hll regions. The FUV and Ha fluxes were corrected for Galac- 
tic extinction using E(B - V") = 0.07 ( [van den Bergh|2000| ) and 
[Cardelli et al.| ( |1989l ) extinction law. The internal Ha extinction for 
each region was obtained using the Ha and 24 /xm fluxes and as- 
suming that the absorbed Ha luminosity in the region scales with 
the 24 /xm luminosity with a factor of a = 0.031 ( Kennicutt et al.| 
|2007| ). FUV extinction was obtained using the empirical relation 
between the ratio of the TIR to FUV luminosities and the UV spec- 
tral index /3 following the formalism given in Calzettr| ( [200T] ). The 
TIR emission was obtained from the linear combination of 8 /xm 
and 24 /xm fluxes given in Table 1 in |Boquien et al.| p010| ) corre- 
sponding to regions with 45" aperture size. We compare the Ha 
and FUV extinctions in the left panel of Fig. [3] For most of the H II 
regions the relation between both extinctions agrees well with the 
expected relation (A(FUV)= 1.8xA(Ha)) from Calzetti's redden- 
ing curve. The values of the FUV extinction are within the range 
obtained by |Verley et'aLl ( [2009] ). 

2.3 H II region sample in the Large Magellanic Cloud 

Our H II region selection for the LMC was determined by the sam- 
ple presented in |Bell et al.|([2002[), which repre sents a subsample 
of isol ated regions from |Caplan & Deharveng| ( [1986[ ). |Bell et aL] 
([2002) presented the fluxes at Ha and 1500 A for a set of H II re- 
gions in the LMC within a diameter aperture of 4.9^ corresponding 
to 65 pc at the distance of LMC (45.9 kpc; |Fitzpatrick et al.|2002] ). 
We performed the morphological classification using the LMC Ha 
image from| Kennicutt & Hodge| ( [T986] ) with an angular resolution 
image of ^7 arcsec. The classification was done visually by one 
of the authors (MR). We applied the same criteria as those for the 
LVL and the M33 sample and catalogued the Hll region sample. 
We found 20, 16 and 16 compact, mixed and shell regions, respec- 
tively. 

Foreground Galactic extinction correction was applied to both 
Ha and FUV (1500 A) fluxes using the Galactic reddening E(B- 
V)=0.06 given in [Bell et al.| ( |2002| . These authors also estimate 
the internal Ha and FUV extinction for each region. The former is 
based on the Balmer decrement and the latter on the TIR/FUV ratio 
with the calibrations derived in|Gordonet al. (2000). We compare 
both extinctions in the right panel of Fig. [3] Most of the regions 
show higher values of FUV extinctions than those derived using 
Calzetti's extinction curve. However, as we study here the H II re- 
gions from M33 and LMC in a separated way the fact that we apply 
different extinction laws for the two samples does not affect the fi- 
nal results of this paper. 



3 RESULTS 



In Fig. [4] we reproduce figure 5 from |Lee et aL] ( |2009| ) using a 
colour-code for each classified galaxy: compact are blue, mixed 
are green and red are shells. In general, the galaxies classified 
as mixed are slightly more luminous than the compact or shells: 
(log L (Ha/ergs-^)) = 38.8 =b 0.2, (log L (Ha/erg s -^)) = 
39.3 =b 0.2, (logL(Ha/ergs"^)) = 38.7 ± 0.2 for compact, 
mixed, and shells, respectively. This is expected as the mixed 
galaxies are formed by several intense knots of star formation in- 
tertwined with some diffuse emission. The mean value for the 
log(SFR(Ha)/SFR(FUV)) ratio in galaxies classified as mixed is 
higher than for the compact and shells, the shells having the lowest 
ratio: (log(SFR(Ha)/SFR(FUV))) is -0.32 ± 0.06, -0.22 ± 0.06 
and —0.43 ±0.07 for compact, mixed and shells, respectively. In an 
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Figure 1. Ha images for examples of compact, mixed and shell objects: galaxies in the LVL sample are shown in the top row, while the star-forming regions 
in M 33 are shown in the bottom row. 
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Figure 2. FUV images of the same objects as those shown in Fig.[T] galaxies in the LVL sample are shown at the top, while the star-forming regions in M 33 
are shown in the bottom row. 



attempt to quantify this trend we have performed a linear fit to the 
data in Fig.|4]for each classification. The fits are shown in the figure 
with the corresponding colour line for each morphology: blue for 
compacts, green for mixed and red for shells. The results of the fit 
are presented in Table 1 . The linear fit for the shell morphology has 
the steepest slope, while the mixed show the shallowest one. We 
will discuss these results in section|4] 

We used the Ha and FUV extinction derived in the previ- 
ous section to correct the observed Ha and FUV fluxes of the 
Hll regions in M33. Then, we converted them into SFRs using 
the relations given in |Kennicuttl ( |1998^ ). These relations assume a 



constant SFRs over ^^lOOMyr and are suitable for galaxies rather 
than for Hll regions caused by short-lived bursts of star forma- 
tion. However, we still apply these relations in order to compare 
the results with the LVL sample. Since we are using the ratio be- 
tween the two SFR calibrators our conclusions will not change if 
we use other calibrations more suitable for H II regions. In Fig. |5] 
we present the logarithm of SFR(Ha)/SFR(FUV) as a function of 
the logarithmic Ha luminosity for the Hll regions in M33. The 
colour code is the same as in Fig. |4] The mixed regions are in 
general more luminous than the compact regions and the shells: 
mean values for the logarithmic Ha luminosity for each class are: 
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Figure 3. Left: FUV extinction versus Ho; extinction for the Hii region sample in M 33. The continuous line represents the relation A(FUV)= 1.8x A(Hq;) 
derived using the Calzetti's reddening curve. The Ha extinction was obtained using the Ho; /24/im ratio. Right: Same figure but for the H ii region sample in 
LMC. The Ha extinction was obtained using the Ha /H/3 emission line ratio. In the lower-right corner of each plot we show the difference between A(FUV) 
and the values expected following Calzetti's reddening curve. 
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Figure 4. Logarithm of the SFR(Hq;)/SFR(FUV) versus logarithm of the 
SFR(Hq;) for the LVL galaxy sample (from Fig. 5 in Lee et al. 2009) with 
colour-code in terms of morphology. Blue: compact galaxies, green: galax- 
ies catalogued as mixed, and red: galaxies with Ho; shell morphology. A 
least-square linear fit is shown for each morphology. The open and com- 
pact black diamonds correspond to the standard and minimal models of the 
IGIMF from ,Pflamm-Altenburg et al.|p009) . 



(log L (Ha/ergs-^)) - 37.6 ± 0.1, (log L(Ha/erg s -^)) = 
38.1 =b 0.2 and (log L (Ha/erg s "^)) = 37.6 ± 0.1 for compact, 
mixed and shells, respectively. The mean values of the SFR ratios, 
(log(SFR(Ha)/SFR(FUV))), are 0.07 ± 0.06, 0.03 =b 0.03 and, 
—0.03 it 0.03 for compact, mixed, and shells, respectively. The 
difference between (log(SFR(Ha)/SFR(FUV))) for compact and 
shells is 0.10 ±0.07, a marginal difference of \Ag. The continuous 
lines in Fig. [5] show the linear fit to the data in each classification 
(see Table 1 for fitting coefficients). 

Ha and FUV fluxes for the regions in the LMC have been cor- 
rected from the extinction derived in ^Bell et al. (2002) and SFR 
calibrations have been applied to derived the corresponding SFRs 
using the relations in 'Kenni cuttj { [19981 ). The results are shown in 
Fig. [6] as well as the linear fits for each classification (see Table 1 
for fitting coefficients). For this sample, there is also a marginal sep- 
aration between the log(SFR(Ha)/SFR(FUV)) ratio for compact 
and for shell regions: (log(SFR(Ha)/SFR(FUV))) are 0.27 ± 0.06, 
0.19 ± 0.08 and 0.15 =b 0.07 for compact, mixed and shells, re- 



Figure 5. Logarithm of SFR(Hq;)/SFR(FUV) versus logarithm of the Ha 
luminosity for the H ii regions in M 33 sample. Extinction corrections to 
Ho; and to FUV luminosities have been applied. A least-square linear fit is 
shown for each morphology classification. 



spectively. The difference in the mean values between compacts 
and shells is 0.12±0.09, a 1.3cr difference. 

We have not combined the H II regions of M 33 and LMC be- 
cause the samples were not created using the same procedure. For 
M 33 we used the Ha image from |Hoopes & WalterboslpOOO] ) and 
identified those isolated H II regions with clear morphology. For 
LMC, we did not create the sample, we rather performed the mor- 
phological classification of a set of H II regions previously iden- 
tified with the purpose of studying the B aimer extinction in each 
region (see |Caplan & Deharveng|1986] ). Therefore, the H II regions 
in LMC are biased to the most luminous ones emitting at Ha and 
H/3. This bias can produce systematic differences in the Ha lumi- 
nosity ranges between M 33 and LMC samples that can affect the 
comparison when the combined sample is studied as a single set. 

We quantify the distribution of the log(SFR(Ha)/SFR(FUV)) 
ratios using histograms for the LVL, M 33 and LMC samples 
(Fig. [7] upper, mid, and bottom panels, respectively). We only show 
the compact and shell classifications to clarify the results. The bins 
in log(SFR(Ha)/SFR(FUV)) are 0.2 dex for the three sets. The sep- 
aration between compacts and shells for the LVL, M 33, and LMC 
samples are 0.11 ± 0.10, 0.10 =b 0.07, and 0.12 ± 0.09. In the 
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Figure 6. Logarithm of SFR(Hq;)/SFR(FUV) versus logarithm of the Ho; 
luminosity for the H ii regions in LMC sample. Extinction corrections de- 
rived from [Bell et al.|p002) have been applied. A least-square linear fit is 
shown for each sample. 



Table 1. Linear fits for the data: log(SFR(Ha)/SFR(FUV))= 
ax log(L(HQ;))+6. 



Sample 


Type 


a 


h 


correl. coeff. 


LVL 


compact 


0.28 ±0.06 


0.3 ±0.1 


0.70 


LVL 


mixed 


0.18 ±0.08 


0.1 ±0.2 


0.50 


LVL 


shells 


0.33 ±0.07 


0.3 ±0.2 


0.73 


M33 


compact 


0.2 ±0.1 


-6 ±5 


0.40 


M33 


mixed 


0.25 ±0.04 


-9 ±2 


0.67 


M33 


shells 


0.21 ±0.06 


-8 ±2 


0.43 


LMC 


compact 


0.1 ±0.2 


-5 ±6 


0.19 


LMC 


mixed 


0.6 ±0.2 


-21 ±8 


0.58 


LMC 


shells 


0.1 ±0.2 


-5 ±8 


0.15 



three samples we find marginal separations, ^ 1.1 — lAa, for the 
compact and shell distributions. 



4 DISCUSSION 

In the previous section we have presented the trends of the 
SFR(Ha)/SFR(FUV) for each data sample: the dwarf galax- 
ies with moderate SFR from the LVL sample and the set of 
Hll regions in the nearby M33 and LMC. Most of the ob- 
jects in the three samples are in the range of low SFR regime, 
-5 < logSFR(Ha)/(M0yr < -1. In each sample the 
mean values of log(SFR(Ha)/SFR(FUV)) are lower for the objects 
catalogued as shells than for the rest of the objects, either mixed or 
compact regions. For the LVL sample we find a difference in the 
log(SFR(Ha)/SFR(FUV)) for compact and shells of 0.11 ± 0.10 
and for M 33 and LMC samples we find differences of 0.10 ± 0.07 
and 0.12 ±0.09. 

If the shell morphology of star-forming regions favours the ex- 
istence of zones of low density gas via holes or fragmentation of the 
shell walls, and the ionising photons are indeed more likely to leak 
the star-forming regions between the fragments, then the difference 
in log(SFR(Ha)/SFR(FUV)) between compact and shells could be 
due to the leakage of ionising photons. 

It is interesting to note that [Hoopes & Walterbos] ( |2000[ ) did 
not find a correlation between the diffuse and ring Hll regions 
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Figure 7. Histograms of the log(SFR(HQ;)/SFR(FUV)) for the three sam- 
ples. All data are extinction corrected. Red histograms are for shells and 
blue for compact regions. 



in M 33 and the fraction of ionising photons escaping the re- 
gions. However, these authors are restricted to regions in M 33 
for which they were able to study the stellar content. These Hll 
regions have low luminosity (only a few are more luminous than 
10^^ erg s~^ after extinction correction has been applied) and they 
represent a small fraction of the H II region population in M 33 (see 
fig. 1 of |Hoopes & Walterbos|200()| . In order to infer firm conclu- 
sions about whether leakage is more important in shells and ring- 
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like regions than in compact objects, a detailed study of the stellar 
content of H II regions with different morphology would be needed. 

Assuming that compact regions are able to absorb all the 
ionising radiation -which would be a strong case, as a compact 
region with a low filling factor would be able to leak radiation 
through the zones between clumps (see [Giammanco et al.||2004| 
|2005| )- and the shells leak all the ionising radiation coming from 
the stellar cluster; then we could say that the difference between 
the log(SFR(Ha)/SFR(FUV)) of compact and shells would be en- 
tirely due to a leakage effect. The maximum observed difference of 
0.12 between the compact and shell regions translates into a fac- 
tor of 0.76 between the SFR(Ha)/SFR(FUV) ratio. Therefore, we 
can say that leakage could change the SFR(Ha)/SFR(FUV) ratio 
between to 25 per cent. 

We can compare our results with those expected from model 
stellar populations. We use the single- star stellar populations ofEl-" 
^dridge & Stanwayi(200 9 ) that assume the H II region stellar popu- 
lation is described by an instantaneous burst with a Chabrier IMF 
( |Chabrier||2003] ). These predict a log(L(Ha)/L(FUV)) of the or- 
der of 0.4 for an H II region 4 Myr old. However, for our observed 
H II regions, log(L(Ha)/L(FUV))?^0.01-0.02, much lower than the 
model ratios. Therefore, leakage of ionising photons is expected if 
the age of the H II regions are ^ 4 Myr. Indeed, for NGC 604 and 
for the set of Hll regions from |Oey & Kennicuttl ( |1997j ) included 
in our LMC sample we observe log(L(Ha)/L(FUV))-0.01-0.02, 
consistent with a leakage fraction (^0-50 per cent) estimated for 
these regions (.Eldridge & Relano 2011 for NGC 604, and |Oey &| 
|Kennicutt|1997| for the LMC H II regions). 

The models only predict ratios of the order of 0.01 at an age 
of 10 Myr. It is unlikely for compact Hll regions emitting at Ha 
to be this old, but shell-like H II regions are expected to have ages 
greater than 4 Myr and the assumption of no leakage could be plau- 
sible. Recently, I Whitmore et al.| ( [20TT] ) have shown evidence of the 
relation between the region morphology and the age of the central 
stellar cluster: very young (a few Myr) clusters show the Ha emis- 
sion of the ionised gas coincident with the cluster stars, for slightly 
older clusters 5 Myr) the gas emission is located in small shell 
structures around the stars, and in still older clusters 5- 10 Myr) 
the Ha emission shows even larger shell structures. They also show 
that if no Ha emission is associated with the cluster then this is 
generally older than ^ 10 Myr. Following this study our sample 
of objects would be younger than 10 Myr, as they still exhibit Ha 
emission associated with the stellar clusters. Moreover, assuming 
a standard shell expansion velocity of 60kms ~^ the radii of the 
shells at t = 10 Myr would be ?^600pc (Relano & Beckman|2005l 
|Verley et al.|2010| ). Most of the observed H II regions in our sample 
are smaller, therefore their age must also be less than lOMyrs. 

We have run models to describe the behaviour of the 
log(SFR(Ha)/SFR(FUV)) ratio as a function of age. For a 10^ M© 
stellar cluster, typical of the most luminous Hll regions, the 
log(SFR(Ha)/SFR(FUV)) ratio decreases by a factor of more than 
one order of magnitude between the young (2-3 Myr) and the more 
evolved (5-8 Myr) H II regions. For a coeval stellar population we 
would then expect a separation between compact and shells of one 
order of magnitude if age is the only responsible to explain the dif- 
ferences in log(SFR(Ha)/SFR(FUV)). 

Other mechanism that can decrease the number of ionising 
photons is the absorption of Lyman continuum photons by the dust 
inside the regions (McKee & Wimams|1997[|Hirashita et al.|2003l 
[Iglesias-Paramo et al.,|2004| ). This mechanism will decrease the 
SFR(Ha)/SFR(FUV) ratio, as less ionising photons will be able to 
ionise the hydrogen. We would then expect that regions with lower 



SFR(Ha)/SFR(FUV) would have higher dust fractions; namely, 
shell regions would have higher fractions of dust than compact re- 
gions. This seems not plausible as shells normally have swept up 
the gas and dust and do not show in general emission at 24 /xm (see 
IVerley et al.|20T0t . 

Despite of the intrinsic difficulties in quantifying the escape 
fraction of ionising photons from star forming regions, we show 
here that this mechanism has an effect in the observed trends for 
SFR(Ha)/SFR(FUV). Leakage mechanism, in combination with 
possible variations of the IMF: vi a stochasticity dCervino & Luridi-| 
janaj2004|, random sampled IMF ([Corbelli et al.|2009 l or assuming 
an IGIMF ( jPflamm-Altenburg et al.|2009| , might well explain the 
differences in the SFR(Ha)/SFR(FUV) for star-forming objects. 



5 SUMMARY AND CONCLUSIONS 

We have studied the SFR for a sample of Hll regions in M33, 
LMC, and a set of dwarf galaxies in the LVL sample. We have 
classified them in compact, mixed, and shell regions and analysed 
the SFR(Ha)/SFR(FUV) ratio for the three morphology types. For 
M 33 and LMC samples we obtained differences of 0.10 ±0.07 and 
0.12 zb 0.09 in the (log(SFR(Ha)/SFR(FUV))) between compact 
and shell regions, a separation of ^ 1.4 — 1.3a. For the LVL sam- 
ple we found a smaller difference, (log(SFR(Ha)/SFR(FUV))) = 
-0.32 zb 0.06 for compacts and (log(SFR(Ha)/SFR(FUV))) = 
—0.43 ± 0.07 for shells, corresponding to a difference of 0.11 ± 
0.10. The maximum observed difference of 0.12 between the com- 
pact and shell regions translates into a factor of 0.76 between the 
SFR(Ha)/SFR(FUV) ratio. 

We present here observational evidence that the escape of ion- 
ising photons from individual star-forming regions in M33, LMC, 
and a sample of dwarf galaxies might account for up to a 25 per 
cent decrease in the SFR(Ha)/SFR(FUV). Rather than performing 
a detailed quantification of the leakage fraction of ionising pho- 
tons we show here that this mechanism should be considered when 
trying to explain the lower than expected SFR(Ha)/SFR(FUV) 
for low star-forming objects. It will probably be a combina- 
tion of different mechanisms: leakage, variations of IMF, SFH, 
and age effects, which might finally explain the decrease of the 
SFR(Ha)/SFR(FUV) ratio for low luminosity systems. 
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